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Abstract. We present a time-dependent Ginzburg-Landau model of nonlinear elasticity 
in solid materials. We assume that the elastic energy density is a periodic function of 
the shear and tetragonal strains owing to the underlying lattice structure. With this new 
ingredient, solving the equations yields formation of dislocation dipoles or slips. In plastic 
flow high-density dislocations emerge at large strains to accumulate and grow into shear 
bands where the strains are localized. In addition to the elastic displacement, we also 
introduce the local free volume m. For very small m the defect structures are metastable 
and long-lived where the dislocations are pinned by the Peierls potential barrier. However, 
if the shear modulus decreases with increasing m, accumulation of m around dislocation 
cores eventually breaks the Peierls potential leading to slow relaxations in the stress and 
the free energy (aging). As another application of our scheme, we also study dislocation 
formation in two-phase alloys (coherency loss) under shear strains, where dislocations glide 
preferentially in the softer regions and are trapped at the interfaces. 
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1. Introduction 

Nonequilibrium states under shear have been extensively studied for near-critical 
fluids and for various complex fluids including polymers, liquid crystals, surfac- 
tant systems, colloidal suspensions and so on [1]. In the soft condensed matter 
physics, considerable attention has also been paid to jamming rheology observed 
in sheared states of supercooled liquids, soft glassy materials such as dense mi- 
croemulsions or granular materials [2]. In these systems while the particle size 
ranges from microscopic to macroscopic lengths, universal constrained dynamics is 
realized under external shear forces. We mention that mesoscopic dynamic het- 
erogeneity and strong shear-thinning behavior have been observed in supercooled 
liquids (at relatively high T) [3-5] and dense microemulsions (at effectively low T) 
[6]. In engineering, on the other hand, plastic flow has long been studied in crys- 
talline and amorphous solids and in glassy polymers [7-10]. In crystals irreversible 
motions of dislocations give rise to plastic deformations and large strains produce 
high-density dislocations. In amorphous solids shear strains tend to be localized in 
narrow shear bands in plastic flow above a yield stress. The width of such shear 
bands is microscopic in the initial stage but can grow to mesoscopic sizes, sometimes 
resulting in fracture. Shear bands were numerically realized at large shear strains 
in molecular dynamics simulations of two-dimensional (2D) two-component glasses 



[11,12] and in simulations of a 2D phenomenological model [13]. Shear bands were 
observed also in granular materials (at effectively zero temperature) [14]. 

We also stress relevance of dislocations in various phase transformations in solids. 
Since the first work by Cahn [15] most theoretical studies have been focused on the 
coherent case in which the lattice planes are continuous through the interfaces 
[1,16,17]. In the incoherent case dislocations appear around the interfaces and the 
continuity is lost partially or even completely. Such incoherent microstructures 
emerge in various alloys when the lattice constants or the crystalline structures of 
the two phases are not close [18,19]. Moreover, they are produced in plastic flow 
because dislocations generated tend to be trapped at the interfaces [20] . Mechanical 
properties of two-phase solids are hence very different from those of one-phase solids. 
Simulation studies in this area are still at the beginning [21-23]. 

In this paper we will present our recent nonlinear elasticity theory to understand 
these diverse effects in 2D [24-26] . Our approach is phenomenological, but its merit 
over microscopic molecular dynamics simulations is that we can put emphasis on 
any aspects of the phenomena by controlling the parameters or changing the model 
itself. In Section 2 we will present a first version of our dynamic model, where 
the gross variables are the elastic displacement u and its velocity ul = du/dt. 
The stress-strain relation and emergence of plastic flow at large strains will be 
discussed. Metastability of the resultant defect structures will be ascribed to the 
Peierls potential [8,27]. In Section 3 we will introduce a free- volume field m [28] 
into our dynamic model, which will drastically change the dynamics of dislocations 
on long time scales. As a result, we can gain insights into the physical mechanism of 
aging in amorphous solids [29-33]. We expect that m is a key order parameter for 
amorphous solids or glasses. In Section 4 we will include the composition difference 
ip = ca — cb for binary alloys to numerically realize dislocation networks in two- 
phase states under shear strain. Such complex effects have rarely been studied in 
physics but are of great technological importance. 

2. Mechanical model 

2.1 Nonlinear elastic energy 

We present a nonlinear elasticity model in 2D [24,25]. In terms of the displacement 
u = (u x ,u y ) from a reference crystal state, we define the strain components as 

ei = X7 x it x + V y u y , 

e 3 = ^xUy + W y u x , (2.1) 

where V x = d/dx and V y = d/dy. We call e\ the dilation strain, e 2 the tetragonal 
strain, and e 3 the shear strain. If we suppose a 2D triangular (or square) lattice, 
the elastic energy should be invariant with respect to the rotations of the reference 
frame by ±mr/3 (or ±nn/2) (n = 1, 2, • • •). Under rotation of the reference frame 
by 0, the shear strains e 2 and e 3 are changed to 

e' 2 = e 2 cos 29 + e 3 sin 29, e' 3 = e 3 cos 29 - e 2 sin 29. (2.2) 
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Figure 1. Stress-strain curves from (2.6). Those at 7 = 1(T 4 and 1(T 3 with- 
out defects at t — exhibit pronounced overshoot. Snapshots of Se-s — e 3 — 7 
are given at 7 = 0.3 and 0.4 (right). No marked overshoot appears in the 
other curve at 7 = 10~ 3 , where high-density dislocations are initially present. 



The elastic energy is written as F — J drf c \ with the density in the form, 

/ e i = ^e? + $( e2 ,e 3 ), (2.3) 

where K is the bulk modulus assumed to be a constant and $ is the shear defor- 
mation energy. For a triangular lattice the simplest form of $ is given by 
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3 — cos7r(\/3e3 — e 2 ) — cos7r(\/3e3 + e 2 ) — cos(27re2) 



(2.4) 



which is invariant with respect to the rotation by 7r/6 in (2.2), is a periodic function 
of e3 with period 2/\/3 for e 2 — (simple shear deformation), and becomes the 
usual form ^(e^ + eg) /2 in linear elasticity for small strains. A characteristic feature 
is that <!> in (2.4) is highly isotropic around its minima in the e 2 -e3 plane. That 
is, around the reference state e 2 = e 3 = 0, $ may be teated as a function of 
e = (e| + e 2 ) 1 / 2 only for e < 1/2. For a square lattice structure we have recently 
used [26] 

$ = 0[1 - cos(27re 2 )] + 0[1 - cos(27re 3 )], (2.5) 

which becomes (/x 2 e| +^ 3 e§)/2 for small strains. One of the principal crystal axes is 
parallel to the x axis in (2.4), while it is along or make an angle of 7r/4 with respect 
to the x axis in (2.5). If we allow the system to rotate as a whole, / e i becomes 
dependent on the rotational strain u> = \7 x u y — \7 y u x . In this sense our model 
is still incomplete and further generalization is needed, particularly in describing 
polycrystalline states or melting phenomena. 
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Figure 2. Displacement u around type C and type CC slips in the most 
favorable orientation in shear strain. The arrows are from the original unde- 
formed position to the displaced position. 



2.2 Dynamics 

We assume that the lattice velocity = du/dt obeys 
d 

p-Q t VL = V- a + + r^VV + V- ct r , (2.6) 
where the stress tensor a= {o^} is given by 

&xx = K<Zl + -q^- , &yy = J<ei - -7^-, &xy = &yx = -q^- (2-7) 

It follows the relation V- o = —SF/Su. The mechanical equilibrium condition 
V- <x = is equivalent to the extremum condition SF/Su — 0. We introduce the 
viscosity rjo and the random stress tensor er R , [34]. We integrated (2.6) in 2D on a 
128 x 128 square lattice with the mesh size equal to the lattice constant a. Use is 
made of (2.4) (where a triangular lattice is supposed), but almost the same results 
follow also using <f> in (2.5). Here /i is assumed to be a constant fiQ. Space and 
time are measured in units of a and tq = (p/^o) 1/ ^ 2ft - The dimensionless viscosity 
is t]q = i] /t ho an d is set equal to 1. Spatial derivatives are appropriately defined 
to yield well-defined microscopic slips. 

Fig.l displays the stress-strain curves at constant shear rate 7 applied for t > 
in units of /io and r _1 . For the curves of 7 = 1CP 3 and 1CT 4 with the pronounced 
peak, we set u = at t = supposing a perfect crystal at t — 0. They approach the 
elastic instability point 7 = \/3/6, where da xy /dj = 0. Then the shear stress drops 
sharply with catastrophic slip formation. See the two snapshots of Se^ = e 3 — 7. 
For the other curve at 7 = 10~ 3 wc initially put high-density dislocations produced 
after cyclic shearing [25]. The stress overshoot is suppressed with increasing the 
initial disorder. In addition, the slip spacing depends on 7 in plastic flow. For 
7 = 1CT 4 the spacing is a few times wider than in the case of 7 = 1CP 3 . 




u 21 - u 20 




Slip length 

Figure 3. Slip energy difference AF B \ ip (£) in the range 20 < £ < 21 for 
7 = —0.04, 0, 0.08, and 0.1 in the model without free volume (left). Displace- 
ment difference 1x21 — 1*20 at 7 = needed for growth by unit length (right). 



2.3 Slips and Peierls potential 



A fundamental flow unit in our problem is a slip composed of a pair of edge dislo- 
cations. As shown in Fig. 2, slips in 2D are divided into those of clockwise (C) type 
and those of counterclockwise (CC) type. A slip along one of the crystal axes is 
metastable without applied strains, where the slip length ia is an integer multiple 
of the lattice constant a. It satisfies the mechanical equilibrium condition and can 
be calculated numerically. For noninteger I there arises an increase of the elastic 
energy, Up(£), called the Peierls potential [27]. Furthermore, under applied strains, 
(63) = 7 an d (e-2) — e, the free energy to create a slip is given by [25] 



\n£ 



slip 



2tt(1 - v) 



T (7 cos 2ip - e sin 2<p)£ + U P (£), 



(2.8) 



in units of fia 2 for the nearly isotropic <E> in (2.4). Here v = (1 — /j,/K)/2, ip is 
the angle of the slip with respect to the x axis, and Up(£) vanishes for integer £. 
In front of the second term, — is for type C and + is for type CC. For simple 
shear deformation with 7 > and e = 0, the most favorable slip orientation with 
the lowest F s u p is <p = for type C and ip = ir/2 for type CC as in Fig. 2. For 
uniaxial stretching with 7 = and e > 0, it is given by <p = — 7r/4 for type C 
and ip = 7r/4 for type CC. In agreement with the last relations shear bands under 
uniaxial deformation have been observed to make angles of ±7r/4 with respect to the 
stretched direction in isotropic amorphous solids [10-13] and granular materials [14]. 
In real crystals slip orientations can be strongly influenced by crystal anisotropy. 

For noninteger £, Up(£) is path-dependent. Here we calculate F s \i p (£) in the range 
20 < £ < 21 for a type C slip parallel to the x axis under shear strain 7. Use is 
made of the linearly extrapolated displacement, 

u e = (1 - a)u 20 + au 2 i, (2.9) 

where a = £ — 20. The it 2 n and it 2 i are numerically obtained for £ = 20 and 21. 
Fig. 3 shows the free energy difference (left) 
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Figure 4. Shear stress (u xy ) (left) and defect energy density fo defined by 
(2.12) (right) vs strain ~((t) for cyclic shear deformation at 7 = ±1CT 3 . The 
shear stress vanishes at the three points A, B, and C (x). 



AF slip = i^upW - F slip (20), 



(2.10) 



which exhibits a maximum at an intermediate i. The maximum approaches 20 as 
7 T 7ci and 21 as 7 J. 7 C 2, where 7 c i ~ 0.1 and 7 C 2 <~ —0.1 here. Thus the slip is 
metastable for 7 C 2 < 7 < 7 c i, but it expands for 7 > 7 c i and shrinks for 7 < 7 C 2. 
Fig. 3 also indicates that about 10 particles move significantly when the right end 
of the slip moves by unit length (right). Against thermal agitations slips can be 
long-lived if the maximum of Ai^ii p is much larger than k^T. 



2.4 Cyclic shear and defect energy density 

Next we apply a cyclic shear deformation, where j(t) = 10~ 3 in the time regions 
nt p < t < (n + l/2)t p and 7(i) = -10~ 3 in the time regions (n + l/2)i p < t < 
(n + l)t p . For the first two cycles, Fig. 4 shows the stress-strain curve (left) and a 
defect energy density ,/b (right). The latter is defined as follows. We divide the 
strain j(t) into an elastic strain 7 e i and a slip strain 7 S — 7 — 701- Roughly speaking, 
the elastic strain outside the slip lines gives rise to the average stress, while the slip 
strain is caused by the jumps of u x across the slips. Regarding the shear stress in 
(2.7) as a function of e<i and e%, we define 7 e i by 

(7x3/(0, 7el) = {(?xy)- (2.11) 

We define the average defect energy density by 

/ D = (/el>-$(0,7el), (2-12) 

where $(e2, 63) is given by (2.4). The elastic energy density stored is then the sum 
of the regular elastic energy density and the defect energy density. For small 7 e i 
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we have 7 e i = (c X y)/^o and $ = poj^/2. In Fig. 4 we notice (i) residual strains 
at vanishing stress at points A,B, and C, (ii) no overshoot in the stress from the 
second cycle, (hi) that (<r X y) < at the ends of each cycle, and (iv) that (<T xy )(t) 
and fo(t) are roughly constant in plastic flow. 

For 7(i) < 0.3 in the first cycle, however, fr>(t) represents the elastic energy due 
to the inhomogeneous fluctuations of the local strains (mainly due to 5e^) with the 
peak height at 0.012 (not shown in the figure). After this initial period, fv(t) is 
in a range of 0.004 — 0.005 reasonably representing the elastic energy due to the 
defects produced in plastic flow. Note that /d(£) depends on 7 and it is almost 
constant around 0.002 in plastic flow with 7 = ±10~ 4 . 

3. Dynamic free-volume model 

3.1 Free volume 

In the presence of point defects or local free volume, there can be a difference be- 
tween the normalized density deviation Sp/p and — e\ even in the linear elasticity. 
Cohen, Flemming, and Gibbs [35] constructed a hydrodynamic description of solids 
including a new variable which they called the vacancy field. Its small deviation 
(= 5p/ p + ei) was predicted to relax diffusively (vacancy diffusion). More generally, 
since vacancies are point defects, let m represent the local free- volume fraction de- 
fined at each lattice point [28] , which can take continuous values. Here we introduce 
m in our nonlinear elasticity scheme [24]. 

As a reference state we suppose a closely packed crystal state with vq being a 
unit cell volume (or area in 2D). We consider the cell number density n cc ii, which 
changes with lattice dilation, and we define m as 



where M is the mass in a unit volume of the reference state (=constant). We assume 
that m is non-negative (m > 0) neglecting excess packing such as interstitials. If 



m = w ("coil - p/M) 



(3.1) 
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Figure 5. Free volume m accumulated at the dislocation cores due to 
m-dependence of the shear modulus. 
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Figure 6. Free energy density (/)(i), shear stress {a X y}(t), and free- volume 
variance */ (8m 2 ) (t) vs A shear flow with 7 = 10~ 3 is applied in the time 
region < t < 400 and is stopped at t = 400. Relaxations become apparent 
with increasing the average free volume (m). 



the density deviation dp is small compared with the average density p, m may be 
approximated as 



m = (to) - (ei + Sp/p). 



(3.2) 



The average free volume (to) = fo((^ccii) — p/M) decreases with lowering T and/or 
increasing p and represents the distance of the system from the closely packed 
perfect crystal. If the total volume of the system is fixed, to is a conserved variable 
and (to) is a constant independent of time. We shall see that (to) is an important 
parameter controlling slow relaxations of the defect structures. 



3.2 Free energy and dynamic equations 



In this work we use the following free energy density, 

/ = Amln(m/ (to)) + amei + ^Ke\ + <&(m, e 2 , e 3 ). 



(3-3) 
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Figure 7. Snapshots of free volume m, shear strain deviation Se3, and free 
energy density / at t = 400, 3200, and 2 x 10 5 for (m) = 0.2, showing growth 
of heterogeneity of m, an increase of plastic strain, and a decrease of /. 



The first term is an entropic term, though its form is not well justified. However, 
we may well replace it by the simpler form A'm 2 /2 without loss of essential results 
of our simulations. The coupling parameter a is assumed to be small since \Sm\ <C 
| <5ei | should hold except close to the dislocation cores. We use the nearly isotropic 
$ in (2.4), but the shear modulus p decreases exponentially with increasing m as 

ju(m) = fio exp(— CSm), (3.4) 

where 5m = m — (m) and C is a positive constant. This relation simply arises from 
the requirement that solids should soften with increasing the free volume. The 
same exponential form was used to describe softening due to point defects [36]. If 
\Sp\ <C p and the strains are small, the mass density deviation Sp obeys 

JUp = -pV • v. (3.5) 

The mass velocity v is governed by 
d 

P-q-v = V- a +T]oV 2 v + V- ct r . (3.6) 
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Including the contribution from the pressure deviation p(SF/6p)u the force density 
on the right hand side is written in terms of F = J drf as 

where use has been made of the identities (S/Su) p = (5/Su) m — \7(S/Sm)u and 
(5/5m) u = —p(5/5p)u- The expressions for follow from (2.7) if Ke\ is replaced 
by Kei + am and the derivatives are taken at fixed to. The lattice velocity du/dt 
is different from the mass velocity v and is of the form [35], 

d -u = v^=v-L(m)( 6 - j P\ +C R , (3.8) 



dt \ Su 



where L(m) is the kinetic coefficient dependent on to and £ R is the random velocity 
related to L(m) via the fluctuation-dissipation theorem. These dynamic equations 
are approximate. However, if the random forces are neglected, the total free en- 
ergy Ft = J dr[f + pv 2 /2] decreases in time without external forces. The system 

becomes steady only in mechanical equilibrium V- <r= 0. 
From (3.5)-(3.8) we derive the equation for to, 



J^m = V • L{m) 



,fS_F 



V — - V- a 



u 



V-Cr- (3.9) 



If the couplings to u are neglected, it follows the nonlinear diffusion equation 
dm/dt = V • D(m)Vm with D(m) = L(m)A/m. Because the diffusion of to 
should be drastically slowed down for small to, we assume the following form, 

D(m) = D exp(-B/m), (3.10) 

where D and B are positive constants. This form is also not well justified. It is 
worth noting that the same form was used for the relaxation rate of the free volume 
in granular matters [37]. Another aspect is the presence of V- o in the equation of 
to. A similar stress-diffusion coupling is well-known in the dynamic equation of the 
composition for gels and viscoelastic fluid mixtures [1,38,39]. See the last section 
for more comments. 



3.3 Slow relaxation due to free-volume accumulation 

We choose v, it, and to as fundamental variables and solve (3. 6), (3. 8), and (3.9) 
by setting A = 0.05^ , B = 1, and D = 2 x 10~ 3 [40]. In Fig.5 we first show a 
steady-state profile of m around a slip for (to) = 0.2, C = 3, and a = 0. Here the 
steady-state condition SF/Sm =const. yields 

mocexp[(C$-aei)/A]. (3.11) 

Accumulation of to around the dislocation cores becomes signficant for small A and 
large C. Such heterogeneity of to should cause slow structural relaxations or aging, 
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Figure 8. Slip energy difference Aj p s ii p (£) in the range 20 < I < 21 at 7 = 
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(left). The quasi-static condition 8F/Sm =const.is assumed along the path 
connecting the steady slip solutions with I = 20 and 21. 



as expected in the literature [29]. As an example, a shear flow at rate 7 = 1CP 3 is 
applied in the time region < t < 400 and is stopped at t = 400. Fig. 6 displays 
time-evolution of the average free energy density (f)(t), the average stress (a xy )(t), 
and the variance \J (8m 2 )(t). The last quantity represents heterogeneity of m. The 
four curves are those for (m) = 0.15, 0.2, 0.25, and 0.3. The other parameter values 
are given in the figure. The relaxations are more apparent with increasing (m). 
We are adding very small Langevin noises [34], so noisy behavior appears only at 
very long times. Therefore, the relaxations here are due to deterministic diffusive 
accumulation of m around randomly and densely distributed dislocations. If the 
noise amplitude is increased, the relaxations become faster and noisy. 

The origin of the slow relaxations may be ascribed to weakening of the Peierls 
barrier with accumulation of m at the dislocation cores. In Fig. 8 the slip energy 
AF s iip defined in (2.10) is calculated in the range 20 < t < 21, where u is fixed 
as in (2.9) but m is relaxed until V5F/Sm nearly vanishes at each a. Namely, F 
is minimized with respect to m at fixed u. The slip-length change supposed here 
is quasi-static taking place only on time scales longer than a 2 /D(m). From Figs. 3 
and 8 we recognize that the maxima of AF s i; p much decrease with increasing (m) 
and that the metastability is broken by applied strains much smaller than in the 
absence of free volume. 



4. Two-phase alloys 



4-1 Composition coupled to nonlinear elasticity 

We consider a binary alloy composed of A and B atoms neglecting point defects. 
The compositions, ca and cb, of the two components satisfy ca + Cb = 1- In 
real metallic alloys an order-disorder phase transition can also occur together with 
phase separation [1,16,21]. However, the composition difference ip — ca — cb is the 
sole order parameter in this work. The free energy density / is of the form, 
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C K 
I = /bw(VO + 2"|V?M 2 + ae^ + —e\ + $(^>, e 2 , e 3 ). (4.1) 



The first term is the Bragg- Williams free energy density [1], 

(4.2) 



JBW 



-^ln(l + ^) + i-^ln(l-^) 



2v 



where vo is the volume of a unit cell, T is the mean-field critical temperature with- 
out coupling to the elastic field. For deep quenching this form is more appropriate 
than the Landau expansion form. The second term in (4.1) is the gradient free 
energy. There are two couplings between \p and the elastic field in /. First, a is 
the strength of the coupling to ei arising from a difference in the atomic sizes of 
the two species, leading to lattice misfit in two-phase states. In particular, around 
a dislocation it gives rise to a compositional Cottrell atmosphere [7] in one-phase 
states and preferential nucleation [42] in the metastable region. Second, the shear 
deformation energy $ is of the form of (2.5) with 

/"fc = Mfco + WaV" (fc = 2,3), (4.3) 

where fj,ko and fiki are constants. If /121 > and /U.31 > 0, regions with larger 
(smaller) ip are harder (softer) than those with smaller (larger) ip. The inhomo- 
geneity in the shear moduli, called elastic inhomogeneity, plays a decisive role in 
the domain morphology in late stages, because it gives rise to asymmetric elastic 
deformations in the two phases [1,43,44]. 

In the coherent case the mechanical equilibrium condition V- 0= may be 
assumed even in dynamics. Here the stress tensor a— {a^} is expressed as in (2.7). 
Using this condition in the linear elasticity, the elastic field has been expressed in 
terms of ip in the previous theories [1,15-17]. We then find the following: 
(i) The typical strain produced by the lattice misfit is given by 

e = aA^/2L , (4.4) 

where Aip is the difference of ip between the two phases and L (= K + (i 2 o in 2D) is 
the longitudinal modulus. In the coherent condition we require eo + 7 < 1/4 under 
applied strain 7. (ii) In the weak limit of cubic elasticity and elastic inhomogeneity, 
one-phase states become linearly unstable for T < T s /(1 + (t/;) 2 /2), where 

T s = T Q + v Q a 2 / L k B (4.5) 

is the spinodal temperature of one-phase states at the critical composition. 

In the incoherent case the mechanical equilibrium does not hold when dislocations 
are created and are moving. We thus set up the dynamic equation (3.6) for u. Then 
the damping of u is governed by the kinetic viscosity r]o /p. The composition obeys 

where u is fixed in SF/5ip. We neglect Langevin noise terms in integrating these 
dynamic equations. In the dilute limit ca — * 0, ca relaxes diffusively with the 
diffusion constant D = XoksTv^ 1 . 
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Figure 9. Strain amplitude e (left) and free energy density / (right) around 
an incoherent hard domain. Here e is zero within the hard domain and non- 
vanishing outside. The peaks of / are located at the dislocation cores. 

4-2 Dislocation formation 

We integrate the dynamic equations in 2D on a 256 x 256 square lattice under the 
periodic boundary condition by setting K//i2o — 4.5, a/^20 = 0-6, AbTo/voA^o = 
C/a 2 /i2o = 0.05, and /U21 = M3i = 0.6/i20 [26]. Space and time are measured in 
units of a and To = (p//i2o) 1 ^ 2 fl, respectively. We set the diffusion constant Do 
at 5 x 10 -5 ?7o/p and the dimensionless viscosity J?o/ T oM20 a t 0.1. The time scale 
of ip is much longer than that of the elastic field by four orders of magnitude. In 
real solid alloys, these two time scales are even more separated, probably except 
for hydrogen- metal systems where the protons diffuse quickly [1] . 

Fig. 9 shows an incoherent hard domain obtained at deep quenching at T/T = 1. 
Slight cubic anisotropy /i3o/M20 = 1-1 is assumed, so the domain shape is a rounded 
square. We displays the strain amplitude e = (e| + e^) 1 / 2 (left) and the free energy 
density / in (4.1) (right). The slips make angles of ±7r/4 with respect to the x axis 
where \e^\ is large, while they are parallel to the x or y axis in the corner regions 
with large 1 63 1 . In the interface region / exhibits a cliff-like structure arising from 
0.08 ■ ■ ■ 1 ■ ■ 1 




Strain Y = Yt 

Figure 10. Stress-strain curves in two-phase states with (ip) = —0.4, 0, and 
0.2. There are no initial dislocations. Yield occurs with dislocation formation. 
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Figure 11. Composition difference tp (left), shear strain deviation 
8e% — ez — 7 (middle), and shear deformation energy $ (right) in plastic 
flow at 7 = 0.2 in Fig. 10. The average order parameter (tp) is -0.4, 0, and 0.2 
from above. The black regions (right) represent harder domains, which are 
enclosed by the softer phase (in gray) . 



the gradient term and higher peaks arising from the dislocation cores. Here we 
remark that dislocation clouds have been observed around 7' precipitates of the 
LI2 structure in alloys with large elastic misfit [19]. 

Next, in the isotropic condition (j,20 — M30 ; we initially prepare coherent two-phase 
states at T = 2Tq < T s for three average compositions. We then apply a shear flow 
at 7 = 10~ 3 for t > 0. As can be seen in the stress-strain curves in Fig. 10, yield 
occurs at 7 ~ 0.1 where dislocations start to be created. In Fig. 11 slips glide in 
the softer regions and stop at the interfaces. We can see that the black points in $ 
representing dislocation cores are mostly trapped at the interfaces. These results 
are consistent with metallurgical experiments [18,20]. 



5. Summary and concluding remarks 

We have shown that the periodicity of the elastic energy density with respect to e 2 
and e3 gives rise to proliferation of dislocations under large strains. Furthermore, 
we have found new approaches to unexplored complex effects with introduction of 
a new variable, the free volume or the composition, coupled to the elastic field. We 
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give some additional remarks. 

(i) In our simulations slips emerge as long straight lines. If disorder is fully in- 
troduced, glide motions of slips in particular directions should be much limited. 
In molecular dynamics (MD) simulations of two component glasses, the degree of 
disorder should be sensitive to the size ratio of the two species. When it was rather 
close to 1, long slips indeed emerged [11]. Performing MD simulations with various 
size ratios is thus informative. 

(ii) We mention the effect of elastic interaction among dislocations. In our case 
slips are more easily created around preexisting ones, leading to growth of shear 
bands, as already reported in Ref. [13]. 

(iii) In Section3, (3.6), (3.8), and (3.9) constitute a two-fluid model. We expect that 
the fluctuations of m and hence Sp can be increased anisotropically under applied 
strains and after migration of m [24]. This effect is intensified for small A in (3.3) 
and large C in (3.4). In polymers and gels the composition fluctuations have been 
observed to increase under strains [1,39]. 

(iv) In Section 4, the composition has been taken as a single order parameter. Ex- 
tension is needed to more general phase ordering involving an order-disorder phase 
transition [1,16,21] or diffusionlcss (Martensitic) structural phase transitions [1,16]. 

(v) We conjecture that the origin of the dynamic heterogeneity observed in MD 
simulations [3-5] could be ascribed to long-range correlations of the free volume. 
This is because the structural relaxation is caused by the free- volume redistribution 
around elastic inhomogeneity. 

(vi) Finally we stress that our nonlinear elasticity is still incomplete. For example, 
there should be the convective terms in (3.9) and (4.6). Also see the comment 
below (2.5). Construction of a more complete theory is under way. 
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